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We propose the pf-\-pd mixing model for CeSb to explain the stability of the ferromagnetic 
double layer structure in the magnetic ordering. The pd mixing causes the saddle type singular 
points, neighboring the A axis, for the bands which gain energy through the pf hybridization 
with the occupied / state. The peak of the density of states due to this combined effect of the 
pf mixing and the pd mixing enhances the stability of the double layer structure. The same 
combined effect also causes the saddle type singular points in the joint density of states of the 
optical transition. The peak structure of the magneto-optical spectra which has been observed 
in experiments is explained by the present model. 
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§1. Introduction 

Cerium monopnictides {CeXp, Xp—P, As, Sb, Bi) are 
semimetals with the NaCl type lattice. Their electronic 
bands have hole surfaces originated from atomic Xp p 
states around F, and electron surfaces originated from Ce 
5c? states around three equivalent X In the tempera- 
ture (T) and the magnetic field (H) diagram,^' '^^ CeXp's 
show various types of magnetic ordered states. They 
consist of the ferromagnetic layers which stack antifer- 
romagnetically. Many types of stacking sequences ap- 
pear in the phase diagram. Among CeATp's, CeSb^-* and 

CeBi^'^-' have the type-IA antiferromagnetic (-|- -I ) 

phase with nearly full (Jz ±5/2) magnetic moment as 
the ground state of weak H case. Here denotes the 
ferromagnetic layer with the magnetic moment parallel 
to the field, and ' — ' denotes the layer anti-parallel to the 
field. 

The pair of ferromagnetic layers (++ or ) is widely 

observed as a unit of the magnetic structure in CeAp's. 
We call it as the 'double layer structure', hereafter. 

The phase diagram of CeSb has been qualitatively re- 
produced by the ANNNI model. However, it is not 
clear how the parameters of the phenomenological model 
are related to those of the microscopic model. 

The origin of the ferromagnetic layer with nearly full 
moment was explained by the pf mixing model^"^^^ pro- 
posed by Kasuya et al. The model also succeeded in 
explaining the origin of the antiferromagnetic stacking 
of the ferromagnetic layers. However, the stability of the 
double layer structure with sufficient energy gain has not 
been obtained. 

The pd mixing term between Ap p bands and Ce d 
bands is neglected in the energy calculation of the orig- 
inal pf mixing model, because the pd mixing effect dis- 
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appears on the A axis where the pf mixing becomes 
maximum. Later, it was pointed out that the pd mixing 
term plays important rolc"'^^"^^-' to explain the Fermi sur- 
faces of CeSb in the ferromagnetic phase^°"^^' and -I- -I — 
phase.25'26) 

In this paper, we employ the pf mixing model includ- 
ing the pd mixing effect, and examine the stability of the 
double layer structure in CeSb. Hereafter, we call the 
model as the pf+pd mixing model. It is shown that the 
p bands, which are pushed up by the occupied / state to 
gain the pf mixing energy, hybridize through the pd mix- 
ing with the d bands around the A axis. This combined 
effect of the pf mixing and the pd mixing causes saddle 
type singular points neighboring the A axis, and sub- 
stantiates the energy gain due to the pf hybridization. 
The double layer structure is stabilized by the combined 
effect. 

We also study the magneto-optical spectra of CeSb. 
Recently, the magneto-optical spectra in various mag- 
netic ordered states have been extensively studied exper- 
imentally by Pittini et al.'^'^^ and by Kimura et a/.^^"^"' 
The individual magnetic ordered state shows its own 
characteristics in the optical conductivity spectra. The 
spectra will give us fine information on the electronic 
states in the magnetic ordered states and work as a good 
probe to study the characteristic band structure of the 
pf-\-pd mixing model. The peak structures observed in 
the experiments are explained by the saddle type singular 
points in the joint density of states, which are caused by 
the combined effect of the pf mixing and the pd mixing. 

In the next section, we introduce the pf-\-pd mixing 
model. In section 3, the stability of the double layer 
structure is examined. In section 4, the optical conduc- 
tivity in the magnetic ordered state is calculated. Fi- 
nally, we summarize the results obtained by the pf-{-pd 
mixing model in section 5. 
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§2. Model 

Let us choose the z axis along the magnetic polariza- 
tion, then the ferromagnetic layer is in the xy plane. We 
denote the layer which has the polarization parallel to 
the axis by the sign '+', and the layer which has the po- 
larization anti-parallel to the axis by the sign '. The 
magnetic ordering is classified by the stacking sequence 
of layers along the axis. 

We employ the following Hamiltonian, 

H ^Hf + Hp + Hd + Hpd + Hj,f, (2.1) 
H} = Y.{Y.esfl{Ri)ft{Ri) 

I t — K^V 

+ Ufffl{ROURi)fl{Ri)U{Ri)}, (2.2) 
ffp = ^ ^ ep(fc)pj(fe)pt(fc), (2.3) 

= E E e4k)dl{k)d,ik), (2.4) 

fe s=a,/3 

^p'' = E E E ^p''('='^'*)ptW^«w 

s=a,fit=K,,v 

+ h.c, (2.5) 

^f/ = EE E ypfik,t,t') 

I t,t'=K.,U 

X ^^■R^fl{Ri)pt,{k) 

+ h.c. (2.6) 

Here, Hf stands for the / states, and Ri denotes site /. 
The quantities fniRi) and fv{Ri) are the annihilation 
operators of the / electrons at Ri with k and v compo- 
nents of the Fg symmetry, respectively.^'^' These states 
have large magnetic moment along the z axis, and are 
considered as the main components of the ferromagnetic 
layers. ^^'^^' k and v are ascribed to the -|- layers and the 
— layers, respectively. £/ is the / level, and Uff is the 
/-/ Coulomb potential. 

Hp stands for the p bands which have hole surfaces 
around F. pt(fc) denotes the annihilation operator of the 
p band electron with symmetry t. There are six sheets 
of p bands. Among them, we consider only k and v 
components of the Fg symmetry which mainly mix with 
K and V components of the / states. The remaining two 
bands of the Fg symmetry (A and /i) are treated as a 
reservoir, and the two bands of the Tj symmetry are 
neglected because they do not have Fermi surface. 

Hd stands for the d bands which have electron surfaces 
around X. ds{k) denotes the annihilation operator of 
the d band electron with the xy symmetry and the spin 
index s. The xy bands mix with the p bands with k and 
v symmetry in the energy region near the Fermi energy, 
Ep. We retain the xy bands in H4, and the other d bands 
{yz and zx) are treated as a reservoir. In the d bands, 
we use the spin index s {a and f3) because the spin-orbit 
interaction is negligible. Hp4 and Hpf stand for the band 
mixing terms. 

In the Hamiltonian, we have retained only the bands 
which are mainly affected by the magnetic ordering char- 



acterized by the ferromagnetic layers in the xy plane. 

Therefore the cubic symmetry of the paramagnetic phase 
is lost at this stage. The A-symmetry axis along the z 
axis is called as the Az axis, hereafter. 

In the phase diagram, there appear layers which are 
denoted by 'O'.'^' They have small sublattice magnetic 
moment and are considered as the layers with the F7 
symmetry. In this paper, we assume that the Fy layer 
is inert for the band structure because the pf mixing is 
weak for the F7 states. 

We introduce Schrieffer- Wolff transformation"^^' to 
avoid treating the / states as the / bands, and to im- 
prove the precision of the numerical process in the band 
energy calculation. With this transformation for the or- 
dered states, the pf mixing term is replaced as the effec- 
tive p-p term as follows: 

N-l ^ 

^p/ = - E E E ^ 

k t,t',t"=K,v n=0 •' 

X Vp*f{k + nAfe, t, t')Vpfik, t, t") 

x^E-*(o^-"^^ 

1=0 

X pl,{k + nAk)pt''{k), (2.7) 

with Ef = and Afc = (0, 0, 27r/(iVa)). N is the 
period of the magnetic ordering, and Ak is the basic 
reciprocal lattice vector corresponding to the period of 
the magnetic ordering, where 2a is the length of the edge 
of the cube of the fee lattice. AA; is obtained by choosing 
primitive translation vectors, one of which represents the 
magnetic ordering and the other two vectors are (2a, 0, 0) 
and (a, a, 0) in the xy plane of the fee lattice. Of course, 
the wave number vector fe -|- nAfe should be contained 
in the first Brillouin zone of the fee lattice. nt{l) is the 
occupation number of t state on the l-th layer in the unit 
cell of the magnetic ordered state. In eq.(2.7), the terms 
which are written as the single site energy of the / states 
are dropped. They are common for k and ly (and also for 
the other states of Fg). This cubic crystal field splitting 
term due to the pf mixing should be included when we 
compare the energy of the phases with various fraction of 
the Fg states. Since the bases of the bands are restricted, 
only the terms which satisfy the condition t = t' = t" in 
eq.(2.7) appear. 

The band parameters are fitted so as to reproduce the 
band structure of LaSb.^' Table I shows the band pa- 
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rameters used in the present calculation. Sa is the atomic 

energy levels, cti, tti and 8\ are the two center integral 
parameters for the nearest neighbor sites. (J2, 712 and (^2 
are the parameters for the next nearest neighbor sites. 
The next nearest neighbor terms are included to repro- 
duce the band dispersion curve precisely on the axis 
where the p/ mixing becomes maximum. The parameter 
values of the pd mixing term and the p/ mixing term are 
the ones used in the previous studies. In later calcu- 
lation, we control the strength of the p/ mixing effect by 
changing — Here, £/ is measured from the Fermi 

energy of the paramagnetic phase. 

The energy of each ordered state is calculated from the 
density of states (DOS) by the equation 

/Bp 
6.E'{DTi{E') + D^{E'))E'. (2.8) 
-00 

Here, is the DOS of the bands which are treated 
explicitly in the Hamiltonian, and £)r is the DOS of 
the reservoir bands. The reservoir is necessary to give 
the right Fermi energy E-^ and the right carrier num- 
ber. As is shown later, of the phases with double 
layer structure has a peak splitting structure across the 
Fermi energy of the paramagnetic phase [E-p ~ eV), 
and this gives energy gain. Without the reservoir, the 
Fermi energy rises up above the peak splitting structure, 
and the energy gain is not obtained. We calculate the 
DOSs (Db and Z)r) by the tetrahedron method. The cal- 
culated energy depends on the number of meshes in the 
fc-space N^. We extrapolate the energy to the Ny, 00 
limit as a function of 1/.^^;^ from the calculated results 
of Nf^ = 20, 30, 40, 50, 60, and estimate the extrapolation 
error AE of the energy. 

§3. Double Layer Structure 

In this section, we calculate the band energy of each 
magnetic ordered state, and show that the pf+pd mixing 
model gives the stabilization of the double layer struc- 
ture. 

Figure 1 shows the energy difference; betwcicn the + -|- 
00 phase and the -1-0 -I- phase as a hmction of s f. The 
energy of the + + 00 phase, which has the double layer 
structure, becomes lower than that of the -l-O-l-O phase in 
the pf+pd mixing model when the / level e j is shallower 
than —1.5 eV. In contrast, the energy difference between 
the two is small in the p f mixing model. The magnitude 
of the energy gain which is obtained by the pf mixing 
model is comparable to that of the previous work.-'^^^ 

We have also calculated the energy difference between 
the -I- + 0000 phase and the +00 -I- 00 phase based on the 
pf+pd mixing model. The former, which has the double 
layer structure, has lower enc;rgy when the pf mixing is 
strong. The -I- -1-000000 phase has also lower energy than 
the +000-1-000 phase in the pf+pd mixing model.33) The 
magnitude of the relative stability of the + + 0000 phase 
and the + + 000000 phase is about 1/3 of that of the 
+ + 00 phase in the range shown in Fig.l. We use the 
value Sf = — 0.7 eV hereafter to analyze the origin of the 
stability of the double layer structure. 

We note that the t = k component hybridizes only 
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Fig. 1. The energy diflference between the + + 00 phase and the 
+0 + phase i7_|_+oo — ^+0+0 a-s a function of the / level Sf. 
□ corresponds to the pf+pd mixing model, and X corresponds 
to the pf mixing model. The error bars indicate twice of the 
extrapolation error AB. 



with the s = a (up spin) component of the d bands in 
the present band scheme, because the other components 
which can weakly hybridize are dropped in eq.(2.1). The 
K bands and the a bands do not hybridize with the t = v 
and s = [3 (down spin) components. Therefore, we call 
the former part as the +-spin bands and the latter as the 
— spin bands. 

Figure 2 shows the density of states (DOS) of the +- 
spin bands near E^. The DOS for the + + 00 phase and 
that for the +0 + phase are shown. The +-spin bands 
are mainly modified by the magnetic ordering, and the 
— spin bands are not modified because fv is not occupied 
in these phases. The p bands, which are pushed up by 

0.4 0.2 




0.4 

DOS (l/eV/Ce) 



Fig. 2. The DOS of the +-spin bands for the + + 00 phase (solid 
line) and that for the +0 + phase (broken line). Horizontal 
lines are Ey for the + + 00 phase (solid line) and that for the 
+0 + phase (broken line). £/ is set to —0.7 eV. 



the strong pf mixing, strongly mixes with the d bands, 
and the combined effect causes the peak structure shown 
in Fig. 2. Detailed discussion of the combined effect will 
be given later. The peak structure in the +0 + phase 
appears above E-p. Therefore, the combined effect does 
not contribute to the band energy gain. In contrast, the 
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peak structure in the ++00 phase appears below Ef and 
above Ep (the peaks marked by a and b). This causes 
the band energy gain. 

Figure 3 shows the band structure of the +-spin bands 
for the + + 00 phase. A part of the p bands is pushed 




W Xx 



Xz' w 



Fig. 3. The band structure of the +-spin bands for the + + 00 
phase. These bands correspond to the DOS in Fig. 2. Xx corre- 
sponds to (1, 0, 0) and X^, corresponds to (0, 0, 1/2). The S' 
axis connects (0, 0, 1/2) and (0, 1/2, 1). The other points and 
axes are the ones of the fee lattice. 



up above Ep on the Az axis due to the pf mixing with 

the occupied / states. The p bands have gaps on the Az 
axis caused by the ordering of the / states. On this axis, 
the pf mixing effect is maximum, while the pd mixing 
disappears due to symmetry. In contrast, the pd mixing 
term recovers its value steeply when the k point moves 
from the Az axis.^^^ Since the p bands, which are pushed 
up by the pf mixing, have energy close to the energy of 
the dxy bands, the band structure is strongly modified 
in the vicinity of the Az axis. This fact is seen in Fig. 3. 
On the axes Ax and E', we can see band splitting effect 
between the p bands and the d bands. 

The p bands, which are pushed up above Ep, have 
small band dispersion along the direction parallel to the 
Az axis. Then the pd hybridized bands have approx- 
imately cylindric equi-energy surfaces around the Az 
axis. This enhances the peak structure of DOS due to 
the van Hove singularity. We note that the peak a in 
Fig. 2 corresponds to the maximum in Fig. 3 just below 
Ep in the Ax axis (the 3rd band from the low energy 
side). 

Figure 4 shows the band structure for the + + 00 phase 
along the (0, 0, 1/4) — (1, 0, 1/4) axis, which is cal- 
culated by the pf+pd mixing model. Figure 5 shows 
one calculated by the pf mixing model. The energies 
of bands in both cases arc identical at (0, 0, 1/4), be- 
cause the pd mixing disappears due to symmetry. When 
the point moves to the (1, 0, 1/4) direction, the band 
splitting due to the pd, mixing is seen in Fig. 4. The max- 
imum marked by a' in Fig. 4 corresponds to the shoulder 
marked by a' in Fig. 2, and the maximum marked by b in 
Fig. 4 corresponds to the peak marked by b in Fig. 2. We 
note that the point a' has maximum type singular point, 
and the point b has the saddle type singular point. 

We note that the DOS of the pf mixing model does 
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Fig. 4. The band structure for the -|- -t- 00 phase along the 
(0, 0, 1/4) — (1, 0, 1/4) axis in the pf+pd mixing model. 
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Fig. 5. The band structure for the -I- -|- 00 phase along the 
(0, 0, 1/4) - (1, 0, 1/4) axis in the pf mixing model. 



not have fine structure, because of the absence of the 
pd mixing term. The gain of the band energy due to 
the combined effect of the pf mixing and the pd mixing 
around the Az axis is important to stabilize the double 
layer structure. 

Let us consider the origin of the relative stability be- 
tween the + + 00 phase and the +0 + phase. The 
latter phase has the period characterized by the wave 
number vector (0, 0, 1). At the point where the main 
band gap is formed, the d band has energy higher than 
the Fermi energy of the paramagnetic phase. Therefore, 
the main band repulsion due to the combined effect oc- 
curs rather higher energy region. On the other hand, the 
+ + 00 phase is characterized by the wave number vector 
(0, 0, 1/2). The gap is formed near the Fermi surface of 
the dxy bands. When we remove the reservoir from the 
present model, Ep locates just above b in Fig. 4, and the 
stability of the +0 + phase is relatively enhanced. 

Figure 6 shows the magnetic phase diagram under the 
field with the parameter £/ = —0.7 eV. As noted pre- 
viously, the crystal field splitting between the Fg states 
and the Fr states, 



EcEF — E] 
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(3.1) 
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magnetic ordered states are given by the following ex- 
pression, 



Fig. 6. The magnetic phase diagram with the / level £/ = —0.7 
eV. The transverse axis is the field H, and the ordinate axis is 
the crystal field splitting Eqef between the F7 layer and the Fg 
layer (Spg — E^j)- 'Para' means the F7 phase. The energy of 
the phases with period-1,2,4 are compared in this figure. 



is a free parameter of the present model. The energy 
difference between the ordered states with different frac- 
tion of the Fg layer depends on the quantity Eqef- At 
the same time, the energy gain from the Zeeman inter- 
action is proportional to the ferromagnetic moment of 
the phase. We have calculated the energy of each phase 
assuming following equation 

E = Sband + (n+ + n-)EcEF - in+ - n-)H. (3.2) 

Here, n+ and n_ are the fraction of the + layers and 
the — layers, respectively. The Ft layers are assumed as 
non-magnetic layers for simplicity. 

When the energy of the Fg layer is enough lower than 
that of the layer (£^cef < —1-4 meV in Fig.6), the 

+ H phase becomes the ground state in the H = 

case. This phase may correspond to the type-IA phase of 
CeSb. The phase turns into the ferromagnetic phase in 
the present calculation, when the field is increased. The 
boundary between the + -\ phase and the ferromag- 
netic phase is about 6 T in the present calculation, and 
is comparable to the experimental magnitude 4 T.^^ In 
the calculation of the Zeeman energy, the 5- value of the 
Fg states is included as •*— qjUbH with gj = 6/7. 

When the energy of the F7 layer (i.e. the paramagnetic 
layer in the present model) is enough lower than that of 
the Fg layer (£^cef > 2 meV in Fig.6), the F7 phase 
(0000) becomes the ground state. This phase turns into 
the ferromagnetic phase by way of the -|- -|- 00 phase as 
the field increases. 

The phase transition with the parameter ii^cEF > 2 
meV will mimic that of CeP, whose ground state is the 
antiferromagnetic F7 phase. In CeP, the phase succes- 
sively turns into the states which have large fraction of 
double Fg layers when the field is increased.^-* 

§4. Optical Conductivity Spectra 

In this section, we calculate optical conductivity spec- 
tra based on the pf+pd mixing model. The optical con- 
ductivity spectra for the ±1 circular light, a±{E), of the 



/ d'fe E E S{Ej{k) - E,{k) - E) 

X I < aj{k)\ds{k + nAk) >< pt{k + nAk)\ai{k) > ^ 
X w{k + nAk) c5r''(s,i) 

+ y d^fe E E S{Ej{k) - Ei{k) - E) 

X I < aj{k)\pi{k + nAk) >< ds{k + nAk)\ai{k) > p 
X w{k + nAk) 4^^(s, t), (4.1) 

Afe = (0,0, 1^). (4.2) 

where aj{k) stands for the diagonahzed band with the 
band index j, and BZ means the integration in the Bril- 
louin zone. iv{k + nAk) stands for the fe dependence of 
the optical transition probability. The term w(fe -|- nAk) 
is proportional to (cosfcx + cos ky)^, since the present 
p ^ d transition is the inter-site process. However, we 
assume the term as constant, because the saddle type 
singular points which correspond to the peak structure 
of the spectra are located neighboring the Az axis where 
kx ^ and ky ^ 0. s and t represent the indices of the 
d, states [a and (3) and the p states (k and u) respec- 
tively. The terms c5.~^''(s,f) and c'!^^{s,t) stand for the 
selection rule of the optical transition. They have 1 for 
c^'^{a,K), c'^^{P,i'), (fr''^{l3,v), cl~^^(a, Ac), and for 
the other combinations of s and t. We note that the op- 
tical transition matrix appears only between and da 
states, or between p^ and djj states in the present simpli- 
fied model. At the same time, the -I— spin bands (which 
consists of Pk. and da) do not hybridize with the — spin 
bands {p^ and djj). By the facts, the optical conductivity 
is separated into the p ^ d components and the d ^ p 
components as expressed in cq.(4.1). 

Figure 7 shows the calculated optical conductivity 
spectra for the ferromagnetic phase. cr-{E) for various 
Ef cases are shown in the figure. When the pf mixing 
is large enough, one of the p bands is pushed up above 
£^F around Xz- Then, the d ^ p transition around Xz 
comes to be possible, and the transition contributes to 
a-{E). The band in the ferromagnetic phase has a sad- 
dle type singular point of the joint density of states at 
Xz, and the peak structure, which corresponds to the 
transition around this point, appears in a-{E). On the 
other hand, ij+ (E) is essentially equivalent to that of the 
paramagnetic phase consistently with the experimental 
results by Kimura et al.^^'^^^ 

The energy of the peak E^eak in Fig. 7 becomes higher 
as the strength of the pf mixing effect, which is propor- 
tional to —l/sf, becomes larger. At the same time, the 
strength of the peak increases as the intensity of the pf 
mixing effect increases, because the fe-space which con- 
tributes to the peak becomes wider. 
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Fig. 7. The optical conductivity spectra u- (E) for the ferromag- 
netic phase. The spectra for Sf = —1.4, —1.2, —1.0, —0.8, -0.7 
eV are shown. The intensity of the pf mixing is proportional to 



Figure 8 shows the relation between the energy of 
the peak Epeak of the ferromagnetic phase and —l/sf. 
The marks '+' in the figure are roughly on the line of 
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Fig. 8. The energy of the peak Epeak of the ferromagnetic phase 
versus the strength of the pf mixing effect (— The peak 
energy is plotted by '+' sign. The broken line represents the 
energy difference at Xz between the p band Ep and the d band 
Ed- 



Ep{Xz) — Ed{Xz), which is the energy difference be- 
tween the p band and the d band at Xz- 

Using the relation between the energy of the peak and 
— l/sf, we can determine — l/ej from the experimen- 
tal results.^*' We obtain the value of the / level as 
Sf = —1.35 eV from Fig. 8, where the peak appears at 
0.4 eV. In later calculation of the optical conductivity, 
we employ this value. In this case, the energy difference 
in Fig.l becomes about half of that of the ef = —0.7 eV 
case, and this case also holds the condition of the rel- 
ative stability of the double layer structure. The value 
of £/ depends on the magnitude of the pf hybridization 
matrix element. "^^^ On the other hand, the energy of the 
peak directly corresponds to the band energy difference 
Ep{Xz) — Ed{Xz)- Therefore, the present choice means 
that the later calculation is carried out based on the band 



parameters which give the band energy difference 0.4 eV 

in the ferromagnetic phase at Xz- When we use this 
readjustment, the calculated spectral change due to the 
change of the magnetic ordering will scarcely depend on 
the detailed choice of the band parameters. 

We note that the transition around Xz becomes a 
maximum type singular point in the pf mixing model. 
Therefore, we do not have peak structure at this energy, 
if the pd mixing is neglected. We also note that the reser- 
voir bands are not modified by magnetic ordering, and 
their optical conductivity spectra are equivalent to those 
of the paramagnetic phase. The low energy peak struc- 
ture which appears in the ordered states will be scarcely 
affected by the contribution from the transitions among 
the reservoir bands. 

Figure 9 shows the calculated optical conductivity 
spectra for the + -\ — phase with £/ = —1.35 eV. Both 
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Fig. 9. The optical conductivity spectra for the -|- -| — phase. 



(7_|_ {E) and cr_ {E) have a peak and a shoulder structure 
which are located at 0.26 eV and at 0.37 eV. The peaks 
mainly consist of d — > p transition, and the shoulders 
mainly consist of p — > d transition. Both the peaks and 
the shoulders have the lower energy than the energy of 
the peak of the ferromagnetic phase. These behaviors 
are consistent with the observation.^*^) 

In the -|- H — phase, the pf mixing pushes up both the 
-I— spin bands and the — spin bands, and the pd band 
repulsion makes saddle type singular points which cor- 
respond to the structure in Fig. 9. In contrast to the 
peak structure of the ferromagnetic phase, the transition 
around Xz does not correspond to the structure. 

Figure 10 shows the contour map of the interband 
transition energy for the + -\ — phase in the fc-space. 
This map has approximately axial symmetry around the 
axis. The peak of a+ (E) and the shoulder of a- {E) 
at 0.26 eV in Fig. 9 correspond to the saddle type singu- 
lar point shown by 'x' in Fig. 10(a). The former corre- 
sponds to the d — > p transition around the point, and 
the latter corresponds to the p d transition around 
the same point. The peak of cr^{E) at 0.37 eV in Fig. 9 
corresponds to the d —>■ p transition around the saddle 
type singular point shown by 'x' in Fig. 10(b), and the 
shoulder of cr+(-E) corresponds to the p ^ d transition 
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Fig. 10. The equi-cncrgy contour map for the optical transition 
in the + H — phase. The transition of the +-spin bands shown in 
(a) contributes mainly to the peak at 0.26 eV of (E) and the 
transition of the — spin bands shown in (b) contributes mainly 
to the peaJc at 0.37 eV of (7-(E). Both peaks correspond to 
the transition from maxima just below Bp to minima just above 
Ef near X on the S' axis in Fig. 10 of ref.30. The equi-energy 
contours of E = 0.3, 0.4, 0.5 eV (a) and those of E = 0.1, 0.2, 
0.3, 0.4, 0.5, 1.0 eV (b) are plotted. The region where the optical 
transition does not occur is closed by the contour named 'edge'. 
X'^ is (0, 0, 1/3), and the A'^ axis connects F and X'^. The 
saddle type singular points are marked by ' x ' in the figures. 



around the same point. 

In Fig.ll, we show the optical conductivity spectra for 
the + H — phase calculated by neglecting the pd mixing 
term (i.e. the pf mixing model). There appears low 
energy transitions, caused partly by the Brillouin zone 
folding effect. This case does not have peak structure. 
That is, the peaks of the optical conductivity spectra 
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Fig. 11. The optical conductivity spectra for the + H — phase of 
the p/ mixing model. 



originate^ from tlic^ saddle type singular points caused by 
the combined effect of the pj mixing and the pd mixing 
effect. 

Figure 12 shows the optical conductivity spectra for 
the +0- phase with ej = -1.35 eV. (y^{E) and 0-{E) 
are identical, and they have a peak at 0.26 eV. The peak 
structure in the +0— phase and the structure of (t+(_E) 
in the + + — phase shown in Fig. 9 have the same energy. 
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Fig. 12. The optical conductivity spectra for the +0— phase. 
IT+ (E) and a- (E) have the same spectrum. 



It is because the — spin bands in the +0— phase are 
identical to those in the + H — phase. The difference is 
on the centroid of the peak, which is shifted to higher 
energy side in the + H — phase than that of the +0— 
phase due to the contribution from the +-spin bands. 
The experimental result of the spectra for this phase has 
not been reported. We expect that the features shown 
in the figure will be observed. 

Figure 13 shows the optical conductivity spectra for 
the + + 00 phase with = —1.35 eV. There appears a 
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Fig. 13. The optical conductivity spectra for the + + 00 phase. 



double peak structure, and both the peaks at 0.35 eV 
and at 0.43 eV in <T-{E) consist of the d ^ p transition. 
(7+ (E) has the structure corresponding to that in (t_ (E) 
with weak intensity. These results are consistent with 
the observation.''"' 

In Fig. 14, we indicate the saddle type singular points 
neighboring the A'^ axis which correspond to the peaks 
of cr-{E) shown in Fig. 13. Similarly to the intcrband 
transition structure of the + H — phase, the dispersion of 
the transition energy is small in the direction parallel to 
the axis around the saddle type singular points, and 
the cylindrical equi-energy region appears in the + + 00 
phase. 
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Fig. 14. The cqui-cncrgy contour map for the optical transition 
in the + + 00 phase. Upper figure (a), for which there appears 
the peak at 0.35 eV, corresponds to the transition from the third 
lowest band to the fourth lowest band shown in Fig. 3, and the 
lower figure (b), for which there appears the peak at 0.43 eV, 
corresponds to that from the second lowest band to the fourth 
lowest band shown in Fig.3. X'^ is (0, 0, 1/4), and the axis 
connects F and X'^ . The saddle type singular points are marked 
by 'x' in the figures. The equi-encrgy contours of i? = 0.3, 0.35, 
0.4, 0.45, 0.5 eV (a) and those ol E = 0.4, 0.45, 0.5 eV (b) are 
plotted. The region where the optical transition does not occur 
is closed by the contour named 'edge'. 



§5. Summary and Discussion 

We have demonstrated that the combined effect of the 
pf mixing and the pd mixing is important to stabihze the 
double layer structure which is widely observed in the 
magnetic phase diagram of CeSb. A part of the p bands 
is pushed up along the Az axis by the pf mixing with 
occupied / state. These p bands strongly mix with 5ri 
bands around the axis through the pd hybridization. 
This causes several bands which have saddle type sin- 
gular points neighboring the axis. The energy gain due 
to the pf mixing is substantiated by the pd hybridized 
bands. The same bands give the saddle type singular 
points of the joint density of states of the optical tran- 
sition in the magnetic ordered states. Characteristics of 
the magneto-optical spectra for various phases of CeSb 
are reproduced by the calculation. 

The energy of the peak structure of the optical con- 
ductivity in the ferromagnetic state corresponds to the 
excitation energy between the pushed up p state and 
the d state at Xz- Therefore the peak structure in the 
ferromagnetic phase of La^Cei-^Sb is expected to shift 
to low energy side as x becomes large, because the pf 
mixing effect decreases. The peak intensity will also de- 
crease. The situation similar to that shown in Fig. 7 will 
be shown. 

In the present paper, we have employed the simplified 
model to carry out numerical calculation with sufficient 
accuracy for the band energy of various magnetic phases. 
The same simplified model accounts the characteristics 
of the magneto-optical spectra of each magnetic phase. 
However, the bands treated as the reservoir will also con- 
tribute to the spectra in various way. For example, the 
transition between the reservoir and the states which are 
modified strongly by the magnetic ordering will be also 
affected. Such transition terms are not considered in the 



present model. Further studies including total bands are 
necessary in quantitative analysis of the magneto-optical 
spectra. 

Recently, Iwasa et al. have proven by the fine X-ray 
diffraction experiment that the 'paramagnetic layer' con- 
sists of the r7 states. The present model ascribing 
the paramagnetic layer as the T-j states will be justified. 
We have treated the layer as the inert state in the band 
calculation. In the strict sense, the r7 states have hy- 
bridization term with the band states. However, it is not 
large for the states near the Fermi energy The essential 
features of the present result will not be modified even 
when we include the hybridization term of the Fy states. 
We expect that the Kondo effect for the paramagnetic 
state^^-* will not be so important. 

We have shown the stability of the double layer struc- 
ture on the basis of the band calculation assTiming the 
long range periodicity. The states which include the dou- 
ble layer structure seems to be more stable than the 
states with isolated layers even when we consider the 
periodicity of eight layers. Therefore, the explanation 
of the stability of the double layer structure on the basis 
of the local electronic structure, in some sense, is de- 
sirable. In the experiments on CeP, we observe phases 
which have very long periodicity with a small number 
of double layer in the sea of the paramagnetic state. 
Such a study on the isolated double layer is retained in 
the future. Nevertheless, we note that, in the short pe- 
riodicity cases, the present band calculation picture will 
be applicable. 
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